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We study fusion reactions of the ^'*Ni+^''^Sn system using the recently introduced density con- 
strained time-dependent Hartree-Fock formalism. In this formalism the fusion barriers are directly 
obtained from TDHF dynamics. In addition, we incorporate the entrance channel alignment of the 
deformed (oblate) ^"^Ni nucleus due to dynamical Coulomb excitation. We discuss the influence of 
particle transfer and other dynamical effects on the fusion cross sections. Calculated cross sections 
are in very good agreement with data and other calculations. 
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I. INTRODUCTION 

With the increasing availability of radioactive ion- 
beams [l| the study of structure and reactions of exotic 
nuclei has become feasible [2]. In particular, detailed 
investigations of the fusion process are crucial for the 
prediction of heavy-element formation and will lead to a 
better understanding of the interplay among the strong. 
Coulomb, and weak interactions as well as the enhanced 
correlations present in these many-body systems. 

Recently, particular attention has been given to the 
^^Ni-|--'^^^Sn system, where a large sub-barrier fusion en- 
hancement was observed [3, IJ For this system fu- 
sion cross sections were measured in the energy range 
142 MeV < Ec.ni. < 195 MeV. In particular, it was 
found that fission is negligible for Ec.m. < 160 MeV 
and therefore the evaporation residue cross sections have 
been taken as fusion cross sections. The enhancement of 
sub-barrier fusion was originally deduced from compari- 
son with a barrier penetration calculation, using a phc- 
nomenological Woods-Saxon interaction potential whose 
parameters were fitted to reproduce the evaporation 
residue cross sections for the ^''Ni-l-^^^Sn system 0, 
Similarly, early coupled-channel calculations, which are 
known to enhance the fusion cross sections by consider- 
ing coupling to various excitation channels and neutron 
transfer, have significantly underestimated the low en- 
ergy fusion cross sections for the ^^Ni+^'^^Sn system Q. 
Subsequently, more sophisticated coupled-channel calcu- 
lations lead to an improvement for the description of the 
lower energy data. Finally, the inclusion of the neutron 
transfer channels with positive Q value in addition to 
inelastic excitations resulted in the best description to 
date 

The theoretical analysis of the fusion data generally 
involves determination of a phenomenological ion-ion po- 
tential such as the Bass model d, 0] , the proximity poten- 
tial 3, [^, [13, fll |, o r potentials obtained via the double- 
folding method[ll, [l^, [13, [lB| • Subsequently, the actual 
fusion cross section is calculated by either using barrier 
penetration mod els [? j_ 3, [i3| , or the coupled-channel 
method [H, [H, [MlSrtiJ- The latter includes vari- 
ous excitations of the target and/or projectile using the 
coupled-channel formalism [2lj , as well as the inclu- 



sion of neutron transfer, and can be consistently applied 
at energies above and below the barrier [l3|. Effectively, 
the inclusion of each additional excitation leads to a mod- 
ification of the original inert core ion-ion potential, result- 
ing in a series of effective barriers. One common physical 
assumption used in many of these calculations is the use 
of the frozen density or the sudden approximation. In 
this approximation the nuclear densities are unchanged 
during the computation of the ion-ion potential as a func- 
tion of the internuclear distance. Furthermore, the ef- 
fects included in channel couplings are usually based on 
the static properties of the participating nuclei, which 
may accurately represent the early stages of the colli- 
sion process, but are expected to change as the two ions 
strongly interact. While these methods provide a use- 
ful and productive means for quantifying multitudinous 
reaction data it is desirable to include dynamical effects 
and make contact with the microscopic theories of nu- 
clear structure and reactions. 

Recently, we have developed a new approach for cal- 
culating heavy-ion interaction potentials which incorpo- 
rates all of the dynamical entrance channel effects in- 
cluded in the time-dependent Hartree-Fock (TDHF) de- 
scription of the collision process |23|. These effects in- 
clude the neck formation, particle exchange, internal ex- 
citations, and deformation effects to all order, as well 
as the effect of nuclear alignment for deformed systems. 
The method is based on the TDHF evolution of the nu- 
clear system coupled with density-constrained Hartree- 
Fock calculations to obtain the ion-ion interaction poten- 
tial. Preliminary calculations for the ^'*Ni-|-^'^^Sn system 
highlighted the importance of dynamical deformation ef- 
fects [2J] . Here we give a completed study of fusion cross 
sections using this formalism. 

In the next section we will summarize some theoretical 
aspects of the density constrained TDHF theory along 
with methods to incorporate dynamical alignment into 
our calculations, as well as the method used the calculate 
cross sections from the resulting barriers. In Section IIIII 
we present interesting aspects of the reaction dynamics 
and compare our results with experiment and other cal- 
culations. 
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II. THEORETICAL METHODS 
A. Density constrained TDHF method 

In this subsection we give a qualitative description of 
the density constrained TDHF method which is used to 
obtain the dynamical barriers for the ^''Ni+^'^^Sn system. 
Further details of the method can be found in Ref . [1^ . 

The density constraint is a novel numerical method 
that was developed in the mid 1980's [1^, and was 
used to provide a microscopic description of the for- 
mation of shape resonances in light systems [1^. In 
this approach the TDHF time-evolution takes place with 
no restrictions. At certain times during the evolution 
the instantaneous density is used to perform a static 
Hartree-Fock minimization while holding the total den- 
sity constrained to be the instantaneous TDHF density. 
In essence, this provides us with the TDHF dynamical 
path in relation to the multi-dimensional static energy 
surface of the combined nuclear system. Since we are 
constraining the total density all moments are simulta- 
neously constrained. The numerical procedure for imple- 
menting this constraint and the method for steering the 
solution to PTDHF(r,t) is discussed in Refs. [2^[2^. The 
convergence property is as good if not better than in the 
traditional constrained Hartree-Fock calculations with a 
constraint on a single collective degree of freedom. 

In Ref. 23] we have shown that the ion-ion interaction 
potential is simply given by 

V{R) ^ EuciR) - Ea, - Ea, , (1) 

where i?DC is the density constrained energy at the in- 
stantaneous separation R{t), while Eai and Ea2 are the 
binding energies of the two nuclei obtained with the 
same effective interaction. We would like to emphasize 
again that this procedure does not affect the TDHF time- 
evolution and contains no free parameters or normaliza- 
tion. In practice, TDHF runs are initialized with en- 
ergies above the Coulomb barrier and in Ref. [2^ we 
have shown that there is no appreciable energy depen- 
dence to the barriers obtained via the density constrained 
TDHF method. A detailed description of our new three- 
dimensional unrestricted TDHF code has recently been 
published in Ref. 27]. For the effective interaction we 
have used the Skyrme SLy5 force [28^ including all of the 
time-odd terms. 



B. Fusion with alignment 

In general, the fusion cross sections depend on the in- 
teraction potential and form factors in the vicinity of 
the Coulomb barrier. These are expected to be mod- 
ified during the collision due to dynamical effects. In 
addition, experiments on sub-barrier fusion have demon- 
strated a strong dependence of the total fusion cross sec- 
tion on nuclear deformation [30| . The dependence on 



nuclear orientation has received particular attention for 
the formation of heavy and superheavy elements [3l| and 
various entrance channel models have been developed to 
predict its role in enhancing or diminishing the proba- 
bility for fusion [H, [s^ . Recently, we have developed a 
new approach for calculatin g th e effect of nuclear align- 
ment for deformed systems [331 ]. In essence, the proce- 
dure for incorporating alignment into the evolution of 
the heavy-ion collision dynamics is done in two separate 
steps: a) A dynamical Coulomb alignment calculation 
to determine the probability that a given nuclear orien- 
tation occurs at the distance i?(io)j where the TDHF 
run is initialized. The alignment generally results from 
multiple E2/E4 Coulomb excitation of the ground state 
rotational band. The distance R{to) is chosen such that 
the nuclei only interact via the Coulomb interaction, b) 
A TDHF calculation, starting at this finite internuclear 
distance R{to), for a fixed initial orientation of the de- 
formed nucleus. Since the experiments are usually done 
with unpolarized beams, in a full quantum mechanical 
calculation one would have to average over discrete quan- 
tum mechanical rotational bands. In the classical limit, 
this corresponds to averaging over orientation angles. A 
general study of taking the classical limit of the rela- 
tive nuclear motion during a heavy-ion collision which 
includes inelastic excitations of one of the heavy ions in 
the entrance channel has been given in Ref. [2§] . 

The heavy-ion interaction potential between two de- 
formed nuclei depends on the distance vector between 
their centers-of-mass, R, and on the relative orientation 
of their intrinsic principal axis systems which may be de- 
scribed in terms of three Euler angles a, /?, 7 per nucleus, 
i.e. in the most general case we have 

V = y(R,a«,/3(i),7^i),a(^),/3(2),7(2)) . (2) 

Explicit expressions for this interaction potential within 
the double- folding method are given in Ref. [15]. The 
expression for V can be simplified if the intrinsic nu- 
clear density distributions are axially symmetric; in this 
case, the potential does not depend on the Euler angles 
.-y(i)^.y(2) which describe rotations about the symmetry 
axes. If we put, for convenience, the distance vector in z- 
direction, R = Rez, the potential between two deformed 
axially symmetric nuclei has the structure 

V = V{R,l3^^\p'^^\Aa) . (3) 

Finally, if one of the nuclei is spherical, e.g. nucleus (1), 
the potential is simply given by 

y = (4) 

where the Euler angle /J*-^-* describes the direction of the 
nuclear symmetry axis relative to the internuclear dis- 
tance vector. 

Details of the dynamic alignment formalism are pre- 
sented in [33| . We give here a brief summary: For a 
given incident energy Ec.m. we carry out a semiclas- 
sical Coulomb excitation calculation of the dominant 
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collective levels of the deformed nucleus. The energy 
levels and _EL-transition matrix elements for ^^Ni are 
taken from experimental data [s^: E2+ — 1.346 MeV, 
E4+ = 2.610 MeV and M{E2,0+ 2+) = -27.0 e fm^ 
(oblate deformation). The Coulomb excitation calcula- 
tion starts at very large internuclear distances (about 
1500 fm) when both nuclei may be presumed to be in 
their respective ground states and stops at the ion-ion 
separation distance R{to) (about 16 fm). The Coulomb 
excitation amplitudes determine the probability distribu- 
tion of initial orientations. In Fig. [T] we show the differ- 
ential alignment probability as a function of the Euler 
angle P used in our calculations. 
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FIG. 1: Dynamic alignment due to Coulomb excitation of 
"^Ni. Shown is the orientation probability as a function of the 
Euler angle /3 in a central collision at internuclear distances 
R = 1500 fm (dashed curve) and at i? = 16 fm (solid curve). 



In the case of one spherical nucleus (^^^Sn) and one 
deformed reaction partner (^''Ni), the total fusion cross 
section is given by an integral over all orientation (Euler) 
angles, with solid angle element dil — 2nsin(3d(3 



f dP 



(5) 



where dP/d^ represents the alignment probability and 
^■(i^c.m., ^^) is the fusion cross section associated with a 
particular alignment. 



Cross section calculation 



For a consistent calculation of fusion cross sections at 
above and below the barrier energies we have adopted 
the commonly used incoming wave boundary condition 
(IWBC) method [Ullil. Once the ion-ion potential for 



a particular orientation, V{R,I3), is calculated the two- 
body Schrodinger equation becomes 



-h^ d^ L{L + l)h^ 



2n dR^ 



2^i?2 



+ V{R,I3)~E 



ijL{R,P) = Q, 



where fi is the reduced mass and L denotes orbital an- 
gular momentum. IWBC assumes that once the mini- 
mum of the potential is reached fusion will occur, con- 
sequently no imaginary part of the potential is needed. 
In practice, the Schrodinger equation is integrated from 
the potential minimum, Rmim where only an incoming 
wave is assumed, to a large asymptotic distance, where 
it is matched to incoming and outgoing Coulomb wave- 
functions to obtain the penetration factor, Pl{E,P), as 
the ratio of the incoming flux at Rmin to the incoming 
Coulomb flux at large distance. The total cross section 
is given by 



a{E,P) = ^y2{2L + l)PL{E,P) , 

Kr 



(6) 



with fco — \f2yLE. For the numerical implementation we 
have followed the procedure for the coupled-channel code 
CCFUL described in Ref. [l^l and exactly reproduced 
their results for an inert-core potential. 



III. 



RESULTS 



We have carried out a number of TDHF calculations 
with accompanying density constraint calculations to 
compute V[Ri P) given by Eq. ([T]). A detailed description 
of our new three-dimensional unrestricted TDHF code 
has recently been published in Ref 
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For the effec- 
tive interaction we have used the Skyrme SLy5 force [2^ 
including all of the time-odd terms. The code was mod- 
ified to self-consistently generate initial states for ^"^Ni 
with different orientations. All of our TDHF calcula- 
tions were done at an initial energy of E'c.m. = 168 MeV 
and separation i?(to) = 16 fm. As we have reported in 
Ref. [23| the potential barriers obtained from the density 
constrained TDHF method are not sensitive to the initial 
energy (above the barrier). We have tested this again by 
running a few orientations at 158 MeV and 176 MeV and 
did not observe any appreciable difference. 



A. Particle exchange 

In a TDHF collision leading to fusion the relative ki- 
netic energy in the entrance channel is entirely converted 
into internal excitations of a single well defined com- 
pound nucleus. In TDHF theory the dissipation of the 
relative kinetic energy into internal excitations is due to 
the collisions of the nucleons with the "walls" of the self- 
consistent mean-field potential. TDHF studies demon- 
strate that the randomization of the single-particle mo- 
tion occurs through repeated exchange of nucleons from 
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one nucleus into the other. Consequently, the equilibra- 
tion of excitations is very slow, and it is sensitive to 
the details of the shape evolution of the composite sys- 
tem. This is in contrast to most classical pictures of nu- 
clear fusion, which generally assume near instantaneous, 
isotropic equilibration. 
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FIG. 2: Time-dependence of the isovector dipole and isoscalar 
quadrupole moments for the head-on collision of ^''Ni-l-^^^Sn 
system at £c.m. = 168 MeV and (3 = 90°. 

Recently, the importance of transfer of neutrons with 
positive Q value in fusion has been emphasized [5l.[36l.[37j. 
In TDHF this effect manifests itself as the excitation 
of the pre-compound collective dipole mode, which is 
likely when ions have significantly different N/Z ratio, 
and is a reflection of dynamical charge equilibration. In 
Fig. [5] we show the time evolution of the isovector dipole 
and isoscalar quadrupole moments for the head-on col- 
lision of ^''Ni+i^ssn system at E'en. = 168 MeV and 
(3 = 90°. The Fourier transform of these oscillations 
show a 9 MeV isovector dipole peak as well as 4 MeV and 
7 MeV isoscalar quadrupole peaks. For a time-dependent 
charge distribution it is possible to calculate the 7- yield 
as a function of gamma energy [26| , . The asymptotic 
yield integrated over a spherical surface is given by 



he 1 



^(2L + l)|a(L,MI' 



(7) 



where the amplitudes a(L, huj) 



a{L, hw) 



in 



i{2L + iy.\ 



L + 1 
L 



1/2 



he 



L+2 



MLihiu) , 



(8) 

are given in terms of the Fourier transform of the mo- 
ments of the density 



MLit) = j d\r'^YLQ{f)p{v,t) 



(9) 
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FIG. 3: Dipole gamma yield for the head-on collision of 

64jg-j_^i32g^ system at Ec.m. = 168 MeV and /3 = 90°. 



Computation of the yield given by Eq. [7] shows that the 
dominant contribution is from the dipole mode, which is 
shown in Fig. [31 
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FIG. 4: (Color online) Change in neutron (blue curve) and 
proton (red curve) numbers of the left fragment (originally 
®*Ni) as a function of the ion-ion separation coordinate, R. 
The dashed lines are there to emphasize the asymptotic val- 
ues. 

While the above analysis provides a picture of particle 
exchange dynamics over time it is also possible to ex- 
amine particle exchange during the initial stages of the 
collision, namely from the well separated nuclei to the 
time when the minimum separation is reached. In Fig. [J] 
we show the neutron and proton transfer to the left frag- 
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ment, which in our case is initially ^^Ni. As we can see 
in Fig. [5] the number of neutrons and protons increase 
for the left fragment as the ions come into contact. As 
expected the transfer of neutrons starts earlier since they 
are dominant in the surface region. For small R values 
the increase is simply due to charge equilibration as the 
two nuclei have a substantial overlap. However, at a sep- 
aration of 12 fm, which corresponds approximately to 
the top of the potential barrier as discussed in the next 
subsection, we have about two neutrons transferred, and 
at 11 fm we have as many as six neutrons transferred. 
It is also interesting to observe that the large transfers 
in the early stages of the collision happen quickly, indi- 
cated by the sudden jumps in neutron number. Natu- 
rally, the right fragment undergoes just the opposite of 
these trends. 



height of 155.8 MeV and location of the barrier peak at 
approximately Rb — 12.1 fm. The difference for smaller 
R values is due to the use of the point Coulomb in- 
teraction in the model calculation, which is unphysical 
when nuclei overlap. As seen in Fig. [5] for higher /3 val- 
ues the barrier is lowered and the barrier peak moves to 
larger R values, with the lowest barrier having a height of 
150.1 MeV and peaking around Rb — 13 fm. The physi- 
cal picture which emerges from these calculations is that 
for all energies above 150.1 MeV, which is the peak of the 
lowest barrier, the fusion cross section will be dominated 
by lower barriers as the cross section above the barrier 
is substantially larger than sub-barrier cross sections. In 
other words the only experimental data point that ap- 
pears to be truly sub-barrier is the lowest energy point 
at 142.6 MeV. 



B. Dynamical potentials and mass 

In this subsection we discuss the barriers obtained via 
the density constrained TDHF method, and other re- 
lated quantities, such as the effective mass and the re- 
duced mass. We have performed calculations of V{R,l3) 
in A/3 = 10° intervals from /? = 0° to /? = 90°. In 
Fig. [5] we show all of these barriers. The lowest bar- 
rier corresponds to (3 = 90° orientation of the symmetry 
axis of the oblate ^''Ni nucleus with respect to the col- 
lision axis. Each subsequent barrier is reduced by 10° 
up to the highest barrier at /3 = 0° orientation. For 
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FIG. 5: (Color online) Potential barriers, V{R,P), ob- 
tained from density constrained TDHF calculations for the 
®*Ni-|-^^^Sn system. Angle (3 indicates different orientations 
of the deformed ®^Ni nucleus in A/3 = 10° intervals. 



the case of /? = 0° orientation the calculated barrier is 
almost exactly the same as the empirical barrier used 
in Ref. [3] without channel couplings, having a barrier 
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FIG. 6; Effective mass calculated from Eq. (|ll|l for the head- 
on collision of ''"'Ni-l-i^^Sn system at Ec.m. = 168 MeV and 
13 = 90°. 

In addition to all of the dynamical effects included in 
the potentials V{R,P) it is also possible to construct ef- 
fects due to dynamical mass. Since in the density con- 
strained TDHF method the potential is obtained from 
the TDHF evolution by essentially extracting the internal 
excitation energy from instantaneous TDHF solutions via 
the density constraint, the energy conservation becomes 

E,.m. = ^M{R, p)R^ + V{R, P) . (10) 

For a particular initial orientation at asymptotic energy 
^'c.m. we obtain the collective velocity R directly from 
the TDHF evolution and the potential V{R,f]) from the 
density constraint calculations. Thus, the effective mass 
is given by 



M{R,(3) 



i?2 



(11) 
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In Fig. [6] we show the effective mass as a function of 
the ion-ion separation coordinate R for energy i?c.m. — 
168 MeV and 13 = 90°. The sohd curve is the effec- 
tive mass obtained from Eq. pT|) scaled with the con- 
stant reduced mass fi for the two ions. The ratio starts 
from unity at large distances and increases as R gets 
smaller and as the ions pass the top of the barrier as 
indicated by the dotted line in Fig. [S] This is due to 
the fact that the relative velocity, R, becomes small as 
the ions begin to overlap substantially, while the numer- 
ator of Eg . [TT] remains non-zero since the energy is above 
the barrier. After the initial slowdown the ions acceler- 
ate once more before reaching the composite system for 
which i? = 0. A similar behavior is observed for all val- 
ues of f3. Traditionally, the effective mass is obtained 
from constrained Hartree-Fock (CHF) or adiabatic time- 
dependent Hartree-Fock (ATDHF) calculations [s^ and 
shows a strikingly similar behavior to the ones obtained 
through our method. 




FIG. 7: (Color online) Reduced mass calculated for the head- 
on collision of '^''Ni-l-^^^Sn system at Ec.m. = 168 MeV and 
P = 90°. 



In addition to the effective mass it also possible to 
calculate the dynamical reduced mass. It is well known 
that the naive formula for reduced mass given by 



MR) 



Ai{R)A2{R) 
'A,{R) + A2iR) 



(12) 



does not correctly represent the dynamical behavior of 
this quantity, monotonically increasing from its asymp- 
totic value to the predictable final value when Ai = A2 = 
{Ai + A2)/2, i.e. a composite has been formed and half 
the mass is in the right half and other half on the left half 
of the numerical box. An alternative way to calculate the 
dynamical reduced mass is given in Ref. [39| . Here the 
dynamical centers and momenta for the left and right 



halves are calculated via 

R, = Tr{fp,)/Tr{pi) 
Pi = Tr{ppi) , 

where the index i — 1,2 denotes left and right halves. 
The left and right masses are then calculated by 



(13) 



dRjdt ' 



leading to the dynamical reduced mass 



M(i?) 



77117712 
7711 + m2 



(14) 



with R = i?i — i?2. This quantity is plotted in Fig.[7]along 
with the naive expression, /io(i?). While we observe in- 
teresting structure in /i(i?) the overall change in the mag- 
nitude around the barrier region is not large enough to 
substantially alter the fusion cross sections. This is also 
shown in Fig. [6] where we have also plotted the ratio of 
the effective mass to the dynamical reduced mass (dashed 
curve). Therefore, in practice we have not included ii{R) 
in our cross section calculations. 



C. Cross sections 

We have calculated the total fusion cross section as a 
function of energy using the alignment averaged fusion 
formula given in Eq. ([S]). The total fusion cross section 
using the potential barriers obtained from density con- 
strained TDHF calculations directly (with constant re- 
duced mass) is shown by the blue curve in Fig. [H Also 
shown is the the latest coupled-channel calculation Q 
including inelastic excitation of ^'^Ni to first 2+ and 3~ 
states and ^^^Sn to the first 2+ state, as well as two- 
neutron transfer (red dashed curve). We have also in- 
cluded the dynamical effective mass by making the well 
known coordinate scale transformation ISSi 



dR^{^^ydR. 



(15) 



As a result of this transformation all of the effects of the 
dynamical effective mass are transferred to the scaled 
potential while the reduced mass /i remains constant at 
its asymptotic value. This is convenient for a number of 
reasons, one being that the we can use our fusion code 
without any modifications. The resulting cross sections 
are shown by the black curve in Fig. [51 As we see the 
effect of the dynamical mass is to raise the cross section 
at higher energies. 

With the exception of the lowest energy data point the 
calculated cross sections are in very good agreement with 
data and the extended coupled-channel calculation. We 
emphasize again that our density constrained TDHF cal- 
culations contain no adjustable parameters or normaliza- 
tion factors. In the region around 160 MeV our calcula- 
tions over predict the experimental cross section while the 



7 



coupled-channel one slightly under predict. The agree- 
ment at lower and higher energies are excellent. The 
question about the lowest energy data point is still an 
open one. It is our understanding that a new experiment 
is planned to measure this cross section 40] . If this cn- 
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FIG. 8: (Color online) Total fusion cross section as a function 
of i?c.m.. Shown are the experimental data (filled circles), the 
latest coupled-channel calculation 0] (red dashed curve) , and 
our density constrained TDHF cross sections with dynamical 
effective mass (black curve) and with constant reduced mass 
(blue curve). 

hancement is confirmed it would be a challenge since all 
of the existing theories underestimate the cross section by 
a few orders of magnitude. We have shown in our bar- 
rier calculations that this is the only point which is truly 
sub-barrier fusion. On the other hand such a large value 
for the cross section is more consistent with an above 
the barrier or very close to the top of the barrier energy. 
This would indicate that there is a mode of ^''Ni nucleus 



which is not well described by the current microscopic 
interactions, like the Skyrme force. One possibility may 
be triaxiality since ^^Ni is experimentally found to be 
gamma soft. 



IV. CONCLUSIONS 

We have performed density constrained TDHF calcu- 
lations of fusion cross sections for the ^''Ni-|--'^^^Sn system. 
Our results agree very well with the measured data and 
with a state of the art coupled-channel calculation de- 
spite having no adjustable parameters. This indicates 
that many of the reaction dynamics are included in the 
TDHF description of the collision. As we investigate fu- 
sion reactions involving neutron rich and deformed nu- 
clei it is apparent that an understanding of the structure 
of these nuclei is crucial to the description of the reac- 
tion dynamics. For these nuclei various effects, such as 
inelastic excitations, particle transfer, and other dynam- 
ical effects lead to substantial modification of the naive 
potential barrier calculations which assume an inert core 
and no dynamics. Consequently, the definition of sub- 
harrier fusion becomes ambiguous since it is difficult to 
determine the barrier a priori. Present calculations indi- 
cate that the lowest energy data point may be the only 
one truly at a sub-barrier energy. If a new measurement 
confirms this data point it would be of great interest since 
none of the theories can reproduce this cross section, thus 
indicating a fundamental property that is not included in 
any of the calculations. 
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